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ABSTRACT 

We show that the position of the central dark matter density peak may be expected to differ from 
the dynamical center of the Galaxy by several hundred parsec. In Eris, a high resolution cosmological 
hydrodynamics simulation of a realistic Milky- Way-analog disk galaxy, this offset is 300 - 400 pc (~ 3 
gravitational softening lengths) after z — 1. In its dissipationless dark-matter-only twin simulation 
ErisDark, as well as in the Via Lactea II and GHalo simulations, the offset remains below one softening 
length for most of its evolution. The growth of the DM offset coincides with a flattening of the central 
DM density profile in Eris inwards of ~ 1 kpc, and the direction from the dynamical center to the 
point of maximum DM density is correlated with the orientation of the stellar bar, suggesting a bar- 
halo interaction as a possible explanation. A dark matter density offset of several hundred parsec 
greatly affects expectations of the dark matter annihilation signals from th e Galactic Cent er. I t may 
also suppor t a dar k matter annihilation interpretation of recent reports by 
Finkbeiner (20121 of highly significant 130 GeV gamma-ray line emission 
parsec projected] away from Sgr A* in the Galactic plane. 



Weniger (2012]) and Su fc 



rom a region 1.5° (~ 200 



1. INTRODUCTION 

Dissipationless (dark matter only) N-body simulations 
predict a nearly universal density profile of dark matter 
(DM) halos, the so-called NFW profile (Navarro et al. 
19971, which features a central 1/r density cusp. If the 



and at present a purely astrophysical explanation of these 
signals cannot be excluded. 

More recently, however, there have been surprising 
reports of a highly statistically significant line-li ke fea- 
ture at ~ 130 GeV in Fermi data from the GC (Brin; 



cooling and condensation of gas insid e these halos ( White 
fe Rees|1978||Fall fc Efstathiou|1980j) is gradual, th en so- 
called adiabatic contraction" dhJlumenthal et al. 111986 



mann et al.||2012 Weniger|2012 Su fc Finkbeiner|20T2" 



Gnedin et al.|2004 ) will pull DM into the central regions, 



thereby increasing the central DM density and further 
steepening the slope of its density profile. It is thus nat- 
ural to expect the maximum of the DM density to occur 
at the dynamical center of a galaxy, the lowest point of 
its gravitational potential. This expectation has made 
the Galactic Center (GC) a preferred target for indirect 
DM detection efforts searching for an annihilation sig- 
nal (iBergstrom et al.||1998| iGondolo k Silk||1999| lAhal 
- al.HiUOtit IBaliz et' al.pi)d| jDanmnger 'fe the| 
3llaboration||2012|). 



Although such a DM annihilation line should be loop- 
suppressed by a factor ~ 10 2 — 10 4 , it has commonly 
been considered a "smoking gun" of DM annihilation, 
due to the difficulty of producing such a high energy 
line with astroph ysical processes (ho wever, see 



Inian et al 



Aharo- 



cbeincr 



roman c 
IceCube co 



As d ata fro m theT ermi Gamma-ray Space Telescope 
( Atwood et al.|2009 ) has been accumulating, the number 
of studies reporting gamma-ray "excesses" or "anoma- 
lies" from the GC, that could be interpreted as a DM 



|2012[ ). In the analysis of |Su fe Fin 
Q2012p , the line's significance is maximized at Galactic 
coordinates (l,b) — (1.5°, 0°), i.e. displaced from Sgr 
A*, the presumed dynamical center of the Galaxy, by 
about 200 projected parsec in the disk plane. This offset 
has been viewed as a strike against a DM annihilation 
interpretation of the line signal. 

On the other hand, the GC is a dynamical ly and ener 
getically complex region dGenzel et al.||2010|). It harbors 



a supermassive black hole (SMBH) flGTieTet al.|[2008 



. 2C 



Gillessen et al.|2009 ) , which may have been active as re 



annihilation signal , has steadily grown (see e.g. Good- 
enough fe Hooperj 12009) |Hooper fc Goodenoughll^Olll 



Hooper & Linden||2011 



cently as 10 Myr ago, if the giant Fermi bubbles (Dobler 



et al. 


2010 


Su et al. 


2010) 



Abazajian & Kaplinghat 2012 
These analyses have searched for the broad gamma-ray 
continuum signal thought to arise from the decay and 



trom a black hole accretion eve nt; it is rich in massive 
young stars ( Bartko et al.|[20 09), which likely formed in 



a sing le star burst event t> ± 2 Myr ago ( Paumard et al 



2006 1 ; it hosts plenty of highly energet ic compact objects 



hadronizati on of DM annihilation products (Bergstrom 
et al. 1998 ) . While this is expected to be the dominant 
DM annihilation signature, it is unfortunately difficult 
to distinguish it from conventional sources (e.g. milli- 
second pulsars, Abazajian 2011, or cosmic-ray interac- 
tions with molecular clouds, Yusef-Zadeh et al. 2012), 

mqk@astro.berkeley.edu 



radia t ing in X-rays an d gamma-rays ( Muno et al.| 2003 
2005 Abazajian 2011); and on ~kpc scales, its gravi 



tational potential is non-axisymmetric due to the pres 
ence of a stellar bar and a boxy bulge (|Blitz fc Spergel 



1991| |Martinez-Valpuesta fc Gerhard||201ip . Given that 
our Galaxy is baryon-d ommated inwards of ^5 — 10 



kpc flKlypin et aL]|2002[ ), one might expect astrophysical 
processes to modify the underlying DM distribution in 



2 



significant ways. In principle these processes may even 
displace the maximum of the DM density away from the 
dynamical center, thus greatly affecting the expected DM 
annihilation signal from the GC. 

In isolated galaxy simulations, resonant interactions 
between the stellar bar and the DM halo have been 



shown to flatten a central cusp into a core flWeinberg & 


Katz|2002[ Athanassoula 


2003jJHolley-Bockclmann et al. 


2005| Weinberg & Katz 


2007|). These results remain 



not see such strong effects flSellwood|2003 Valcnzucla & 


Klypin||2003 


Colin et al. 


2006f. Regarding an oft-center 



DM density peak, these interactions are interesting be- 
cause they may also induce a "dark bar" a nd other non- 



axisymmetric perturbations in the DM (Athanassoula 



2002] |Ceverino fc Klypin||2007[ ) 



External gravitational perturbations, for example dur- 
ing a merger or a near passage of a satellite galaxy, could 
displace the baryonic center from the D M. Such off sets 



have been measured in galaxy clusters (Allen 1998), in 
which separations between the center of the X-ray emit- 
ting gas and the gravitational center determined from 
strong l ensing can be as large as ~ 30 kpc for relaxed 
clusters flShan et al.|2010|). A study of THINGS galaxies 



(|Shan et al.|2010|) A 

by Trachternach et al.| ( |2008 ) found that offsets between 
the photometric and dynamical centers were less than 
one radio beam width (~ 10" or 150 — 700 pc) for 13 
out of 15 galaxies with well-constrained photometric cen- 
ters. However, two galaxies in their sample, NGC 3627 
(a barred Sb galaxy showing signs of a recent interaction) 
and NGC 6946 (a barred Scd galaxy), exhibit moderate 
offsets between one and two beam widths. 

Supernova- or AGN-driven gas outflows may rapidly 
and non-adiabatically alter the potential in the central 
regions of proto-galaxies, prior to the formation of the 
bulk of their stars. Repeated episodes of such impulsive 
outflows, followed by slow adiabatic re-accretion of gas, 
may irreversibly transfer energy to the DM, flattening 



the central cusp in the proc ess (Read & Gilmore 2005 - 
Pontzen fc Governato|[20L2"|). Although this effect may 



not by itselt produce an off-center DM density peak, the 
resulting cored density profile will be more susceptible to 
perturbations. 

Lastly, the presence of a SMBH has bee n argued to 
lead to the formation of a steep cusp of DM ( Gondolo & 



Silk 1999 Gncdin & Primack 



centeredon the SMBH. This 



2004) in the inner parsec 
would obviously preclude 



any significant offset between the peak in DM annihila- 
tion signal and Sgr A*. 

In this work we report on the search for the presence of 
an offset between the dynamical cen ter and the maximum 
DM density in the Eris simulation ( |Guedes et al.|[20lT ), 
one of the highest resolution and most realistic cosmo- 
logical simulations of the formation of a Milky- Way-like 
barred spiral galaxy. We find evidence for such an off- 
set, at a scale of 300 — 400 pc, c onsistent with the offset 
seen by Su & Finkbeiner (2012). Since such an offset is 
only seen in the dissipational hydrodynamic simulations, 
and not in our collisionlcss purc-DM runs, we suggest 
that baryonic physics is in some way responsible for its 
formation. We examine a number of different physical 
mechanisms, but the limited resolution of this study does 
not allow us to conclusively settle on a single preferred 
explanation. At this stage, we wish to draw attention to 



the possibility of the maximum DM density (and hence 
annihilation luminosity) not being coincident with the 
dynamical center of our Galaxy, commonly associated 
with Sgr A*. We hope that our results will stimulate 
future work, examining other high resolution hydrody- 
namic galaxy formation simulations, and investigating in 
more detail the physical mechanisms that can give rise 
to an offset DM density peak. 

The remainder of this paper is organized as follows. In 
Section 2 we describe the numerical simulations that we 
have analyzed. In Section 3 we present the evidence for 
a DM offset in the Eris simulation. In Section 4 we go 
over several possible formation mechanisms and confront 
each of them with data from the simulations. In Section 
5 we discuss implications for indirect detection searches 
towards the GC, and finally in Section 6 we present our 
conclusions. 

2. SIMULATIONS 

Our analysis makes use of four cosmological zoom-in 
simulations of the formation and evolution of a Milky- 
Way-analog galaxy: one DM+hy drodynamics simula- 
tion, "Eris" ( |Guedes et~aT1|2011| ; and three pure-DM 
dissipationless simulations, ErisDark" (Pillepich et al., 
in prep.), "Via Lactea II" (Diemand et al. 2008), and 
"GHalo" ( |Stadel et al.||2009p . Details ot these simula- 
tions are presented in the listed references, and we only 
briefly summarize their salient features here. Eris and its 
DM-only twin ErisDark were run wit h the N-body+SPH 
code Gasoline ( |Wadsley et al. 20041, but only Eris uti- 
lized the SPH hydrodynamics. Eris resolves the forma- 
tion and evolution of a galaxy with 1.3 x 10 7 high res- 
olution DM particles of mass 9.8 x 10 4 M and a simi- 
lar number of SPH gas particles with an initial mass of 
1.2.1 x 10 4 M Q . Gravitational interactions are softened 
using a cubic spline density kernel of length e so f t = 124 
proper per] Throughout the course of the simulation, 
star particles are created from dense gas (> 5 atoms 
cm ) according to a heuristic star formation recipe with 
a 10% star-formation efficiency. These star particles 
affect surrounding gas through a supernova blastwave 
feedback prescription that injects thermal energy, mass, 
and metals. The simulation results in a realistic looking 
barred late-type spiral disk galaxy, that matches many 
observational constraints on the structure of the Milky 
Way. For example, it has a low bulge-to-disk ratio of 
0.35, falls on the Tully-Fisher relation, has a stellar-to- 
total mass ratio of 0.04, and a star formation rate of 
1.1 Mq yr^ 1 . Eris is the most realistic such simulation 
available today. 

ErisDark is a DM-only twin to Eris, meaning that it 
was initialized with the same phases of the Gaussian ran- 
dom density field, but all of the matter is treated as 
DM, while in Eris 17% is baryonic. ErisDark thus has 
a slightly higher DM particle mass of 1.2 x 10 5 M , but 
employs the same e so ft . A detailed comparison of the two 
simulations will be presented in Pillepich et al. (2012, in 
preparation). Lastly, we also compare results with two 
of the highest resolution pur e-DM simulations eve r per- 
formed, Via Lactea II (VL2, Diemand et al. 2008) with 
a particle mass of m p — 4, 100 M^ and a force s oftcn- 
ing of e 8oft = 40 pc and GHalo ( |Stadel et aL]|2009[ ) with 



1 Forces become Newtonian at 2 e s . 
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TABLE 1 

Four Cosmological Zoom-in Simulations. 



Name 


Code 


m DM [Mq] 


Csoft [PC] 


iV vir M vir [M©] 


i?vir [kpc] 


(a, r- 2 [kpc], p-2 [M kpc 3 ]) 


Eris 


Gasoline 


9.8 x 10 4 


124 


1.9 x 10 7 7.9 x 10 11 


239 


(0.0984, 18.8, 3.84 x 10 6 ) 


ErisDark 


Gasoline 


1.2 x 10 5 


124 


7.6 x 10 6 9.1 x 10 11 


247 


(0.173, 42.5, 8.73 x 10 5 ) 


Via Lactea II 


PKDGRAV2 


4.1 x 10 3 


40 


4.1 x 10 8 1.7 x 10 12 


309 


(0.144, 53.0, 7.78 x 10 5 ) 


GHalo 


PKDGRAV2 


1.0 x 10 3 


60 


1.1 x 10 9 1.1 x 10 12 


267 


(0.137, 61.1, 3.53 x 10 5 ) 



Note. — "Idm is the mass of high-resolution DM particles, e so ft their gravitational softening. N v i T and M v i r are the 
total number and mass of all particles (including gas and star particles in Eris) within the virial radius R v i r , defined as the 
radius enclosing a density of 92.5 times the critical density. The last column lists the parameters of the best-fit Einasto profile 
(ln(p/p_2) = — 2/a [(r/r_2) ct — 1]) for the mean enclosed density (p)(<r) = M(< r)/(47r/3 r 3 ) (see Fig. [I}. 
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Fig. 1. — Radial profiles of the average enclosed dark matter 
density (p)(<r) centered on the location of the potential minimum, 
starting at the radius enclosing at least 20 particles. The vertical 
dotted lines indicate the softening lengths, e so f t = 124, 124, 40, 
and 60 pc for Eris, ErisDark, VL2, and GHalo, respectively. The 
location of the maximum enclosed density (for r > e so ft) is denoted 
with a filled circle. Dotted lines show the best-fit Einasto profiles 
(fitted to 1 kpc < r < i? v i r ), whose parameters are given in Table[l] 

m p = 1, 000 M Q and e soft = 60 pc. Both VL2 and GHalo 
were run with the purely collisionless N-body code PKD- 
GRAV2. The relevant parameters of all simulations are 
summarized in Tabled] 

3. AN OFF-CENTER DM DENSITY PEAK 

We define the dynamical center of the halo to be the 
location of the minimum of the total gravitational po- 
tential, and select for our analysis all particles within a 
radius of 2.5 kpc from the center. At z = there are 
104,781 DM, 13,792 gas, and 4,597,762 star particles in 
this region in Eris, contributing 1.03 x 10 10 , 2.07 x 10 8 , 
and 1.97 x 10 10 M Q in mass. In ErisDark, VL2, and 
GHalo there are 42,611, 1,718,223, and 4,041,566 DM 
particles in this region, for a mass of 5.06 x 10 9 , 7.04 x 10 9 , 
and 4.12 x 1O 9 M , respectively. 

3.1. Enclosed density profiles 

Fig. [T] shows radial profiles of the averaged enclosed 
DM density, (p(< r)) = M(< r)/(47r/3r 3 ), for the four 
simulations. The densities are higher in VL2 than in 
GHalo and ErisDark owing to the larger virial mass and 
concentration of the VL2 host halo, but in Eris even 



higher DM densities are reached at r > 500 pc due to the 
baryonic contraction. Compared to ErisDark the density 
profile is slightly steeper in Eris, and this is also reflected 
in a much higher halo concentration. At even smaller 
radii, baryonic physics produces a roughly constant mean 
density core {2 The radii at which the enclosed density is 
maximized are denoted with filled circles in Fig. [TJ and it 
is clear that in Eris this maximum is significantly offset 
from the dynamical center, while in the dissipationless 
simulations it occurs at close to one e so ft- 

3.2. Three-dimensional localization of the density peak 

To study the DM density offset in more detail, we de- 
posit DM particles using a cloud-in-cell algorithm onto 
a three-dimensional grid with cell width equal to 10 pc 
for Eris and ErisDark (500 3 grid) and 5 pc for VL2 and 
GHalo (1000 3 grid). The grid is oriented in space such 
that the z-axis is aligned with the disk normal in Eris, 
and with the major axis of the prolate density ellipsoid 
for the other cases. We smooth the discretized 3D density 
fields by applying a Gaussian smoothing kernel of width 
a. Fig. [2] shows slices through this grid with a = e so [ t . 
The position of the cell with the maximum DM density 
is indicated with a red cross, and again it is clear that 
it is significantly displaced from the dynamical center is 
Eris. The maximum occurs at (x, y, z) = (260, 30, —90) 
pc, corresponding to an offset distance D a s — 280 pc, 
which is about 2.3 e so ft. In the dissipationless simula- 
tions on the other hand, D Q ff = 78, 7.1, and 7.1 pc, all 
well within one e so ft. 

We have repeated the determination of D Q s for differ- 
ent choices of a, and the results are shown in the top 
panel of Fig. [3j The values of D ff for a <C e so ft are not 
meaningful, since shot noise can lead to unphysical den- 
sity spikes on very small scales. But for a > e so f t the 
conclusions are robust: the DM offset remains below one 
e so ft, i.e. is consistent with no offset, for all three DM- 
only simulations regardless of smoothing, while for Eris 
the offset is greater than e so f t up to a ~ 260 pc. That 
-D ff drops to zero for even higher a is not surprising, 
since any density distribution will appear centered when 
smoothed on sufficiently large scales. 

In the bottom panel of Fig. [3] we show the maximum 
value of the density, p max , for our four simulations. Of 
course p max is much more sensitive to a than D g, since 
larger smoothing includes contributions from surround- 
ing lower density material. At a — e so ft the density peaks 

2 These density profiles are centered on the potential minimum. 
Centering on the point of maximum DM density slightly changes 
the central slope. 
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Fig. 2. — Slices of the central 2 kpc X 2 kpc DM density field 
at z = through the location of maximum DM density. From 
left to right: x-y, x-z, and y-z planes; from top to bottom: Eris, 
ErisDark, VL2, and GHalo simulations. The position of the slice in 
the perpendicular direction is given in the top left corner of every 
panel. The density field has been smoothed with a Gaussian kernel 
of width cr equal to one gravitational softening length e so f t = 124, 
124, 40, and 60 pc, respectively. The z-axis coincides with the disk 
normal in Eris, and with the major axis of the prolate DM density 
ellipsoid in the other cases. The slice thickness is lOpc for Eris and 
ErisDark and 5pc for VL2 and GHalo. The images are centered on 
the minimum of the total potential (marked with a white cross), 
while the location of the maximum density is indicated with a red 
'x'. In the three dissipationless simulations (ErisDark, VL2, and 
GHalo) the offset between the maximum DM density and the total 
potential minimum is less than e ao f t , but it is 2.3 e so f t in Eris. 

a value of 0.84 M pc -3 in Eris, about 25% higher than 
the value (0.67 M Q pc -3 ) at its dynamical center. In Eris- 
Dark /5 max = 0.40 M Q pc -3 , less than half of the peak 
density in Eris. We caution, that the resolution of our 
numerical simulations is not sufficient to resolve the in- 
ternal density structure of the offset peak. With higher 
resolution run the peak-to-center density contrast may 
well be higher. 

3.3. Evolution of the DM offset 

We have performed the same analysis described in the 
previous section on all 400 outputs of the Eris and Eris- 
Dark simulations. Fig. [4] shows that the offset measured 
in Eris at z = is no fluKe, but persists over cosmological 
time scales. At very early times (z > 2) there is no offset 
in either Eris or ErisDark. Starting at z « 1.5, however, 
the DM density maximum in Eris starts to depart from 
the dynamical center. Over a period of about two Gyr 
the DM offset grows to D Q ff « 340 pc (almost 3 e so ft), 
where it remains for the remainder of the simulation. In 
contrast, in the ErisDark simulation D h remains below 
1 Csoft for almost its entire evolution, albeit with occa- 
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Fig. 3. — The dependence of D a g and p m ax on the width a of the 
Gaussian density smoothing kernel. The gravitational softening 
scale e so f t = 124 pc for Eris and ErisDark is indicated with the 
dotted line. In the dissipationless simulations (ErisDark, VL2, and 
GHalo) the DM offset remains well below e so f t for all a. In Eris 
the offset is larger than e so f t out to a = 260 pc (about 2 e so ft). 
The maximum density is about a factor of two higher in Eris than 
in ErisDark, and depends only mildly on cr for cr > e so f t . p ma x is 
larger in VL2 due to its higher halo mass. 
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Fig. 4. — The time dependence of the DM offset from the total 
potential minimum in Eris (blue) and ErisDark (red). The gravi- 
tational softening e so f t (z) = 124 pc is indicated with a dashed line. 
In ErisDark the DM offset remains around or below 1 e so f t for al- 
most the entire simulation, while in Eris it begins to significantly 
exceed e so f t around z = 1.5 and remains at Ri 3e so ft afterwards. 
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sional spikes up to ~ 200 pc. 

In Eris, D Q g fluctuates around its time average (over 
the last 4 Gyr) of 340 pc with a root mean square 
(rms) dispersion of 51 pc. The closest the peak comes 
to the center over this time is D g — 180 pc. The 
peak preferentially lies near the disk plane; its mean 
vertical (perpendicular to the disk plane) displacement 
is only = 64 pc, with an rms dispersion of 46 

pc. The maximum density varies around a mean value 



of (p n 



= 0.84 M Q pc with an rms dispersion of 



0.02 M pc 3 , and reaches minimum and maximum val- 
ues of 0.79 and 0.92 M Q pc~ 3 . 

Outputs in the Eris simulation are spaced ~ 35 Myr 
apart, which is too long to resolve the dynamics of the 
offset peak, given the local dynamical time of ~ 15 
Myr. Nevertheless it is already clear from looking at this 
coarsely time-sampled data that the peak locations are 
not randomly distributed throughout the central region. 
We defer further discussion of the temporal evolution of 
the DM offset and implications for its physical nature to 
Section |U 

3.4. Correlation with DM density profile flattening 

The growth of the DM density offset appears to be 
well correlated with a flattening of the central DM den- 
sity profile in Eris, as demonstrated in Fig. [5| The top 
panels show mean enclosed density profiles in the inner 
region (r < 5 proper kpc) at several output times from 
z = 3 to z = 0. We see two notable differences be- 
tween Eris (left panel) and ErisDark (right). First, the 
DM densities (plotted in proper units) tend to be higher 
in Eris than in ErisDark, which is indicative of "adia- 
batic contraction" . Secondly, while the enclosed density 
in ErisDark continues to increase towards smaller radii, 
the profiles flatten out in Eris, indicating the formation 
of a core. Note that we use the term "core" loosely, indi- 
cating a substantial flattening of the density profile, but 
not necessarily implying a constant density. 

Inner density profiles are notoriously difficult to prop- 
erly resolve in N-body simulations, requiring high force 
resolution, large particle counts, and accu rate time in- 
tegration with sufficient l y small timesteps flPower et al. 
"20031 JZemp et al | [20071 |Diemand et~aLl|20(J8] |Dubinski 
et al.||2009[ ). In Pillepich et al. (in preparation) we in- 
vestigate the influence of baryons on the full shape of 
the density profiles in more detail, including convergence 
studies. We find a convergence radius for the z = 
density profile in ErisDark of 0.9 kpc, corresponding to 
~ 5000 enclosed particles. The square symbols in the top 
panel of Fig. [5] indicate r 50 oo (the radius enclosing 5000 
DM particles) for the different outputs. Note that r^ooo 
is significantly farther out than 3 e so ft (indicated with the 
vertical bars), a commonly used "rule of thumb" for how 
far in density profiles can be trusted. For Eris, rsooo is 
probably overly conservative, since the potential is domi- 
nated by baryons and many more than 5000 gas and star 
particles are found inside of 7-5000- As an intermediate 
case we also mark with circles riooo? the radius enclosing 
1000 DM particles^ 

In the bottom panel of Fig. IS] we show the evolution 
of the logarithmic slope din p/alnr, measured at these 

3 At z = 0, there are over a million star particles within rsooo 
and more than 400,000 within riooo- 
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Fig. 5. — Evolution of the inner DM density profiles. Top: Den- 
sity profiles at z = 3, 2, 1.5, 1, 0.5, and for Eris (left) and 
ErisDark (right). All quantities are plotted in proper units. The 
black dotted lines indicate the best-fitting modified Burkert pro- 
file (see Eq. The location of rsooo, the radius enclosing 5000 
DM particles, is indicated with squares, rsooo corresponds to the 
density profile convergence radius in ErisDark at z = 0. For Eris 
the density profile can be trusted to smaller radii due to its larger 
particle counts, and so we additionally mark riooo with a circle. 
The vertical bars indicate r = 3e ao f t . Bottom: The evolution of 
the logarithmic slope (dlnp/dlnr) measured at different radii: the 
symbols and thick lines correspond to the slopes at rsooo, and for 
Eris we also show the slope at riooo- All lines are boxcar averages 
over 10 outputs. This plot shows that baryonic physics in Eris 
leads to the flattening of the central density profile. This flatten- 
ing appears to be correlated with the growth of the DM offset (cf. 
Fig. Hi. We emphasize that these slopes are not asymptotic slopes, 
and a value of less than —1 does not imply a strongly cusped profile 
all the way to the center, but merely indicates the local slope at 
r 5000i which may lie outside of the scale radius at early times. 

radii. We obtained these slopes by first fitting the en- 
closed density profiles inwards of 5 kpc to a modified 
Burkert profile, 



(P)(r) = 



Po r' c 



r * ( r 2 +r 2)(/3-«)/2' 



(1) 



and then evaluating its logarithmic slope at the radii of 
interest, 

dln(p) 



din 7 



1 + (r c /ry 



This functional form describes the central enclosed den- 
sity profile in Eris and ErisDark much better than a 
constant power law or Einasto profile. Note, however, 
that this fitting function is designed to track the enclosed 



G 



density profile all the way down to e so ft, so way beyond 
the convergence radius. We wish to emphasize that we 
merely use this fit to evaluate the logarithmic slope at 
radii outside of 3 e so ft, and do not attach any physical 
meaning to the asymptotic inner slopes a preferred by 
the fits. 

Fig. [5] clearly shows that Eris develops a cored density 
profile, in the sense that it has a much shallower loga- 
rithmic slope at rsooo than ErisDark. Measured at rsooo 
the log slope in Eris becomes progressively shallower un- 
til it stabilizes at a value of ~ —0.3 at z ~ 1. Note that 
r 5000 is larger than D g(z) at all redshifts. Measured at 
riooo the profile even flattens out completely with a log 
slope of ~ 0. In ErisDark, on the other hand, the slope at 
r 5000 remains close to the NFW value of —1 after z = 1.5. 
Comparison with Fig. [4] shows that the onset and time 
scale of the core formation in Eris is remarkably well cor- 
related with the growth of the DM density offset. This is 
a strong hint that whatever baryonic process drives the 
core formation may also be responsible for the offset DM 
density peak. 

4. NATURE AND ORIGIN OF THE OFFSET PEAK 

The existence of a well defined offset to the peak of 
the DM distribution in Eris is unexpected, and its phys- 
ical nature not immediately apparent. Some possibilities 
that we have considered are: (i) statistical fluctuations, 
(ii) an incompletely dissolved subhalo core, and (iii) ex- 
citation of a DM density wave by the stellar bar or other 
perturber. In the following sections we address these 
possibilities in turn. 

As we have seen, the output cadence of the original Eris 
simulation is not sufficient to follow the dynamics of the 
peak. To help with understanding the physical nature 
of the offset, much higher cadence outputs are desirable. 
To this end we have restarted the Eris simulation from 
output number 393 [z = 0.0151) with a ~ 20 times finer 
temporal resolution (At = 1.43 Myr) and evolved the 
simulation for 140 time steps (206 Myr). Visualizations 
of nine consecutive high cadence outputs, spanning 11.5 
Myr, are shown in Fig.pJ and a movie of the full evolution 
can be viewed at [http://vimeo.com/45114776. 

Visual inspection of these much higher cadence outputs 
shows that the offset peak appears to sometimes jump 
around somewhat discontinuously between two subse- 
quent outputs. Occasionally there are multiple peaks of 
roughly equal density. As already seen with the coarser 
time resolution, the point of maximum density appears 
to avoid the very center and remains close to the disk 
plane. 

4.1. Statistical fluctuations ? 

As shown in Sec. |3.3[ the growth of the DM offset co- 
incides in time with the formation of a core in the Eris 
DM distribution. In a constant density core, sampled 
with a small number of N-body particles, Poisson fluctu- 
ations can give rise to a spurious density peak offset from 
the center. The mean DM density in the Eris' core (in- 
wards of 1 kpc) is 0.51 M Q pc -3 , corresponding to an N- 
body particle density of 5,300 particles kpc -3 . At z = 
there are 509 particles within 2 e so ft of the location of 
maximum density, much higher than the expected 306 
for the average core density and well beyond what Pois- 
sonian fluctuations can produce. We have numerically 
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Fig. 6. — 3D contour plots of the DM density in the central 
(2 kpc) 3 volume of Eris. The outer contour (light blue) corre- 
sponds to p = 0.45 Mq pc~ 3 , the middle contour (light orange) 
to 0.8 Mq pc - 3 , and the opaque innermost contour (red) to the 
99 th percentile of the DM density in the volume. The dynam- 
ical center is marked with a black dot, and the stellar disk lies 
in the X-Y plane. The images show the time evolution over 11.5 
Myr (At = 1.43 Myr) starting from z = 0.0151. In the last two 
frames we have plotted the positions of all 78 particles located 
within 1 £ so ft of the maximum density at the second-to-last output 
shown. Only 1.43 Myr later (last frame) they have already dis- 
persed throughout the plotted volume, indicating that the density 
offset is not a bound structure. 

confirmed this by constructing 500 randomly drawn par- 
ticle distributions with the same mean particle density 
as the Eris core, and analyzed these samples in the same 
manner as described above. 

In Fig. [7]we compare the distributions of p max and D g 
between the random samples and the last ~ 124 Eris 
outputs spanning 4 Gyr of evolution. Eris' /9 max distri- 
bution is peaked at 0.85 M Q pc -3 , way beyond the ran- 
dom distribution, which peaks at 0.61 Mq pc -3 and has 
an rms dispersion of only 0.01 M Q pc -3 . Similarly, the 
-D ff probability distribution function for Eris is peaked 
around 350 pc and looks nothing like that of the random 
samples, which slowly increases towards larger -Doff, in 

1/3 

accordance with a simple volumetric D oS scaling. These 
comparisons very clearly demonstrate that the DM off- 
sets we have found in Eris are not due to statistical fluc- 
tuations. 

4.2. Incompletely disrupted subhalo core? 

In this potential explanation, the offset density peak 
would be identified with the tightly bound central re- 
gions of an incompletely disrupted subhalo. The subhalo 
would have to have been massive enough to have experi- 
enced substantial baryonic condensation and associated 
contraction of its DM, making it more resilient to com- 
plete tidal disruption than its DM-only counterpart in 
ErisDark. If such a halo fell into the host at z > 1.5, 
dynamical friction would quickly drag it in to the cen- 
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Fig. 7. — Probability distribution functions of the maximum den- 
sity pmax (top panel) and its distance to the potential minimum 
D a ff (bottom), compared between the Eris simulation (blue) and 
random samples (orange). The Eris sample consists of the last 
~ 124 outputs spanning 4 Gyr of evolution. The random sample 
is comprised of 500 realizations of a randomly drawn particle dis- 
tributions with the same uniform mean density as the Eris core 
inwards of 1 kpc. 

ter, where tidal interactions would have stripped most of 
the weakly bound material in its outskirt. The inspiral 
would have stalled once its mass dropped to the point 
where dynamical friction becomes unimportant. After 
that point it would continue to orbit around the dynam- 
ical center, its temporal cohesion possibly aided by the 
stabilizing effect of a h armonic potential (Kleyna et al. 
2003] |Read eTaT|[2006| . 

We can rule out this possibility, based on several obser- 
vations: (i) as mentioned above there is no stellar coun- 
terpart to the DM offset, and (ii) the peak persists for 
hundreds of dynamical times, which seems too long for 
an orbiting subhalo core to survive. The final nail in the 
coffin (iii) is the fact that the offset DM peak does not 
appear to be a bound feature. This is demonstrated in 
the last two panels of Fig. [6j in which we show the posi- 
tions of all 78 particles located within f e so f t of the max- 
imum density peak at t = 10 Myr (after z — 0.015). A 
mere 1.4 Myr later, these same particles have dispersed 
throughout the central volume. This implies that the 
offset density peak is comprised of different particles at 
different times. 

4.3. External perturber? 

Next we consider an external perturber as a possible 
mechanism to give rise to the DM offset. The stellar disk 
in the center is self-bound, in the sense that it does not 
require the DM halo to be held together gravitationally. 
An external perturber, for example in the form of a pass- 
ing satellite, could then impart a kick to the stellar disk 
that might cause it to slosh around in the underlying sta- 
tionary DM halo. In this situation the DM density peak 
would correspond to the cuspy center of the DM halo, 




I 



Fig. 8. — Visualizations of the motion of the potential minimum 
in Eris in absolute (code) units, color coded by redshift. Top: From 
z = 1.5 to z = 1.1. At z = 1.3 a satellite passes near to the center, 
but it does not noticeably perturb the dynamical center. Bottom: 
the last 40 outputs, from z = 0.1 down to z = 0. The location 
of the maximum DM density is marked with a small black sphere 
connected with a thin line to the position of the potential minimum 
at that time. 



and the offset would simply reflect the displacement of 
the potential minimum (the stars). One might expect dy- 
namical friction to quickly damp out any such sloshing, 
but the efficiency of dynamical friction is strongly re- 
duced in a constant density (harmonic) core ( Read et al 
2006 1 , and thus the center of the self-bound stellar disk 



may be able to "orbit" around the DM density maximum 
for many dynamical times. 

Some fraction of the DM (the low velocity tail) would 
also be bound to the disk and would move with the stars. 
In this picture, the existence of a DM offset may be very 
sensitive to the stellar mass of the galaxy. For a DM halo 
of a given mass, a higher stellar mass galaxy would bind 
more of the DM to it, and reduce the strength of the DM 
offset or eliminated it altogether. Hydrodynamic galaxy 
formation simulations suffering from a strong baryonic 
overcooling problem may thus not be able to observe DM 
offsets. 

In Eris we have indeed identified a satellite passing 
close to the center (D ss 65 kpc) around z = 1.3, which 
may be a good candidate for an external perturber, since 
the DM offset begins to grow around that time. The 
satellite's mass is 1.8 x 10 10 M Q at infall (z = 2.7) and 
2.8 x 10 9 M Q at the time of close passage. We have ex- 
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Fig. 9. — Probability distribution functions of the cosine of the 
angle between the direction towards the DM offset and the ori- 
entation of the stellar bar, defined as the major axis of the stellar 
density ellipsoid inwards of 1 kpc. The three panels show the distri- 
butions split by cosmic time: early times, prior to the appearance 
of the DM offset (top panel), an intermediate time period over 
which the DM offset grows in magnitude (middle), and late times, 
during which the DM offset remains more or less constant at ~ 340 
pc (bottom). The early distribution is consistent with uniform in 
|cos 9\, whereas once the DM offset becomes pronounced the distri- 
butions prefer large values of |cos 9\, indicating alignment between 
the offset peak and the stellar bar. 

amincd the motion of the potential minimum in absolute 
code coordinates (top panel of Fig.[8|, and do not see any 
evidence for a sharp kink or abrupt displacement at the 
time of the satellite's passage. Furthermore, it appears 
that the potential minimum is moving smoothly in code 
coordinates, and it is the DM offset that moves around 
the potential minimum, not the other way around (bot- 
tom panel of Fig. [8]). 

4.4. Density wave excitation by the stellar bar? 

Another possibility is that the DM offset is the result of 
a density wave excitation of some kind. In this explana- 
tion the peak would be contributed by different particles 
at different times, just as observed, its spatial evolution 
would be set by whatever external source is providing the 
excitation, and it may be long lived, provided a contin- 
uous excitation mechanism exists. A natural candidate 
for such an external source would be the stellar bar in 
Eris. The potential is completely dominated by the stars 
in the central region, and so it seems plausible that the 
DM distribution may be affected by departures from ax- 
isymmctry in the stellar component. Indeed, resonant 
interactions between the DM halo and a stellar bar have 
been invoked to explain the transformation of an initially 
cuspy DM densit y profile into a cored one ( Weinberg & 
|Katz||2002 " 



2007) 



In isolated galaxy simulations it has been possible 
to identi fy DM halo particles trapped in different res 
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Fig. fO. — Top: Frequency scan of the angular evolution of the 
direction towards the DM offset. R 2 is the sum of the squared 
differences between the best-fit sinusoid of frequency uj and the 
angular position of the DM offset. Lower values of R 2 indicate a 
better fit, and the deep spike to R 2 < 1 at ~ f5Gyr _1 corresponds 
to a periodicity of 69 Myr in the angular evolution. Bottom: The 
best-fitting sinusoid (red) overplotted on angular position of the 
DM offset in the last 22 outputs. 



may l ead to the formation of a "dark bar" ( Colin et al i 
Ceverino fc Klypinj 2007 1 , which may be showing 



2006 



up m Eris as a DM offset. The ring-like morphology of 
particles trapped in the co-rotatio n resonance (see e.g. 
Fig. If of Ceverino & Klypin 2007) is reminiscent of the 
distribution of DM in the plane of Eris (see Fig. [2]) . A 
detailed investigation of the resonant structure of DM 
particle orbits in Eris is beyond the scope of this paper, 
but will be pursued in future work. 

Arguments in favor of the bar-driven explanation are 
that the growth of the offset coincides with the formation 
of a DM core, that the offset peak preferentially lies close 
to the stellar disk plane, and that the direction from the 
dynamical center towards the offset density peak appears 
to be correlated with the orientation of the stellar bar. 
This last point is demonstrated in Fig. |9j which shows 
distributions of |cos 6\, the angle between the offset peak 
and the major axis of the stellar density ellipsoid inwards 
of 1 kpc (i.e. the orientation of the stellar bar). 

We have split the distributions up by time epoch: dur- 
ing early times (< 4 Gyr), prior to the appearance of 
a significant DM density offset, the distribution is con- 
sistent with uniform in |cos 9\, corresponding to a com- 
pletely random placement of the density peak. Between 
4 and 8 Gyr the DM offset first begins to become signif- 
icantly larger than e so ft, and the |cos 9\ PDF is strongly 
peaked towards unity. Its mean value of (|cos 6\) — 0.83 
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corresponds to an angle of 34°. At later times, once the 
DM offset has become established, the correlation weak- 
ens slightly, with (|cos 9\) dropping to 0.77, but it still 
remains preferentially aligned with the bar. 

The top panel of Fig. |10| shows the results of a fre- 
quency scan from uj/2-k — 0.1 to 10 3 Gyr" 1 of the angular 
evolution of the direction towards the offset. For every 
frequency ui, we performed a least squares fit to a sinu- 
soid for all coarsely sampled outputs after 8 Gyr. The 
plot of R 2 = J2i[®ott(ti) — 7r/2cos(w(t i — t )] 2 exhibits a 
distinct low R 2 spike at a frequency of uj/2tt = 15 Gyr -1 , 
corresponding to a period of 69 Myr, as well as at several 
of its harmonics (30, 45, 60 Gyr" 1 , etc.). These spikes in- 
dicate the presence of a periodic signal in the offset angle 
evolution, which is consistent with the alignment of the 
offset with the orientation of the stellar bar. In the lower 
panel of Fig.[l0]we show the best-fit sinusoid overplotted 
on the offset angle evolution for the 22 outputs since 13 
Gyr. 

4.5. Numerical Resolution 

At D a g « 3e so f t , the scale of the offset is worrisomely 
close to the force resolution of the simulation. A detailed 
study of the internal density structure of the offset peak 
must certainly await much higher resolution hydrody- 
namic simulations. Nevertheless, the fact that no offset 
is observed in ErisDark, which has the same force resolu- 
tion as Eris, gives us some confidence that the existence 
of an offset DM density peak is no numerical artifact, and 
indicates that baryonic physics appears to be responsible 
in some way. 

We would like to be able to check whether the offset 
remains at the same physical scale (~ 340 pc) even in a 
simulation with a smaller gravitational softening length. 
At the moment we do not have access to such a sim- 
ulation^] but we can go in the opposite direction and 
look in a lower resolution run. We have run ErisLores 
and ErisDarkLores at eight times poorer mass resolu- 
tion (7.8 x 1O 5 M and 9.5 x 10 5 M Q , respectively) and 
with a four times larger gravitational softening length, 
e so ft=495 pc. We have analyzed these simulations in the 
same way as described above, except that all length scales 
have been increased by a factor of four to account for the 
larger e so f t . We CIC-deposit particles onto a grid with 
cell width of 40 pc and smooth the density field with a 
Gaussian kernel with a = e so f t = 495 pc. 

As with the higher resolution simulations, we find that 
the point of maximum density is displaced from the min- 
imum of the potential in ErisLores: Z? ff = 430 pc at 
z = 0.15 and 650 pc at z — 0. And again the DM-only 
counterpart ErisDarkLores has a much smaller offset, 
D off =60 pc. Although D g is somewhat larger in Eris- 
Lores than in Eris in absolute terms, it is clear that the 
offset is not scaling linearly with e so ft : while -D ff ~ 3e so ft 
in Eris, it is only 1-1.5 e so ft in ErisLores. 

We have also checked time steps and energy conser- 
vation in our simulations, as insufficiently short time 
steps can lead to the formation of a spurio us density 
core (IZemp et al. 1120071 IDubinski et al.||2009|). The Eris 



\J Esoft /ft) whe re a is the acc e leratio n of a particle, as rec- 
ommended by Power et al. (2003). This results in time 



steps as small as 0.16 Myr at low redshift, which is three 
times smaller than the fixed time step that Dubinski et al. 
( 2009 1 employed in a simulation with comparable e so f t 
(their model mlOOK) and which they found to be short 
enough to avoid the formation of an artificial core. 

These checks give us further confidence that the DM 
offset we have observed in Eris cannot be attributed 
solely to insufficient numerical resolution or timestep- 
ping. On the other hand, we cannot claim to have 
fully resolved D Q g either - higher resolution studies are 
needed. 

4.6. The role of the central DM core 

It is clear that the central DM core plays an important 
role for the formation of an off-center DM density peak. 
Without a flattened central density profile it may be dif- 
ficult to excite an offset peak, since any non-central den- 
sity enhancement would be overwhelmed by the steeply 
rising central cusp. Indeed, as we have seen above, the 
growth of the DM offset temporally coincides with the 
formation of such a core. We have not yet determined 
what physical processes led to the formation of this core 
in Eris. 

A core formation mechanism that has recently been 
suggested is the repeated removal of large amounts of gas 
from the central regio ns through the action of violent su- 



pernovae explos ions ([Read & Gilmore||2005[ |Pontzen & 
Governato|2012D. Ifth 



rese baryonic outflows abruptly al- 
ter the potential, the DM may respond by flowing out of 
the center, transforming the cusp into a core in the pro- 
cess. Perhaps this DM redistribution could also create 
an offcenter DM density peak. Note, however, that the 
central potential in the Eris galaxy is already dominated 
by its stellar component by the time the core formation 
begins at z ~ 2. Even though Eris does experience a 
star burst triggered by a merger at z ~ 1, the associated 
supernovae feedback is therefore unlikely to substantially 
alter the central potential. Furthermore, this is a single 
star formation ev ent, while the mechanism of Po ntzen fc] 



Governato (12012) requires multiple application of violent 
outflows followed by gradual re-accretion of gas. 

To further check the role of SN feedback on the forma- 
tion of a DM density core and peak offset, we have ana- 
lyzed the last output (z = 0.7) of the ErisLT run, which is 
identical to Eris except that it uses a lower star formation 
threshold (0.1 cm -3 ). This lower SF threshold reduces 
the efficiency of the supernova-driven gas blowout and 
results in a disk galaxy with a larger bulge-to-di sk ratio 



simulation suite uses adaptive time steps that scale as 

4 We are in the process of re-running the Eris simulation with 
e S oft=50 pc, but this is an expensive simulation that will not finish 
for several months. 



and a mor e compact stellar configuration (see Guedes 
et al.|20iT ). The DM distribution, however, appears not 
to have been affected as dramatically. The DM density 
profile in ErisLT is almost identical to the one in Eris at 
z = 0.7, both exhibiting a core with radius ~ 1 kpc, with 
-Doff=350 pc in ErisLT and Z3 o ff=360 pc in Eris. 

The absence of a dependence in the DM distribution on 
the star formation threshold suggests that SN-driven out- 
flows are not the only core formation mechanism, high- 
lighting a possible important difference between the evo- 
lution of Milky- Way s ized galaxies and that of the d warf 
galaxies described by Governato et al. (2010 2012 1. In- 
stead, perhaps the stellar bar may be responsible for 
this transformation, through resonant angular momen- 
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Fig. 11. — DM annihilation luminosity "surface density" 
(J Pq M dl) in the central 2 kpc X 2 kpc region of Eris. The con- 
trast between the dynamical center (white plus) and the peak (red 
cross) is ~ 11%. The vertical axis corresponds to the disk nor- 
mal, and the coordinate system has been rotated to maximize the 
angular offset of the peak. 

tu rn exchange between the bar and the DM as suggested 
by iWeinberg & Katzl (12002 " 
fed" 



20071. Clearly, a more de- 
he DM core formation process in 



tailed investigation of 

Eris is needed, and will be the topic of a future study. 

5. DISCUSSION 

Provided the central DM distribution in the Eris sim- 
ulation is representative of that in our Milky Way, the 
results we have presented here have several important 
consequences for indirect DM detection efforts aimed at 
the GC. 

The first concerns the shape of the DM density pro- 
file used to predict the annihilation luminosity from the 
central regions of the Galaxy. Outside of ~ 1 kpc the 
DM density profile in both Eris and ErisDark (as well as 
VL2 and GHalo) is well described by an NFW or Einasto 
profile. However, cooling and condensation of gas in Eris 
has dragged DM towards the center, such that the nor- 
malization of the density profile in Eris is larger than 
in ErisDark: at 1 kpc the mean enclosed DM density is 
2.7 times higher in Eris (0.52 M pc~ 3 ) than in ErisDark 
(0.20 Mq pc -3 ). At angles greater than ~ 7 degrees from 
the GC, one may thus expect baryonic processes to lead 
to a significant enhancement of the surface brightness of 
diffuse radiation arisi ng fr om DM annihilation. 

As discussed in Sec. |3.4| at even smaller radii baryonic 
physics leads to the formation of a DM core in Eris, and 
as a result the mean enclosed DM density in Eris never 
rises to more than 0.73 M pc~ 3 . For comparison, an 
extrapolation of the best-fitting Einasto profile for Eris- 
Dark reaches 1 M pc" 3 at 200 pc, 10 M pc~ 3 at 3 
pc, and asymptotes to 90 M pc -3 . We thus caution 
against extrapolating density profiles determined from 
dissipationless DM-only simulations (such as Via Lactea 
II, GHalo, and Aquarius) all the way in to the GC in 
order to infer central annihilation luminosities. 

Secondly, a maximum in the DM density that is not 
coincident with the dynamical center of the Galaxy sig- 
nificantly modifies expectations for where the central sur- 



face brightness of DM annihilation radiation should peak. 
The time-averaged value of D D g in Eris is 340 pc, which, 
if seen edge on, would correspond to an angular offset 
of 2.4 degrees. For comparison, the search for a DM 
annihilation signal by the H.E.S.S. imaging Atmospheric 
Cherenkov Telescope was restricted to a region of 45 —150 



pc in projected distance from the GC ( Abramowski et al. 
20111. In this context, the recent report of a highly sig- 



nificant detection of a ga mma-ray line at 127 ± 2 GeV in 
Fermi data from th e GC ( Su &: Finkbe iner 2012) is par- 
ticularly intriguing. Su fc~Finkbcincr| ( |2012) report that 
the signal is maximized at a Galactic longitude of £ ~ 1.5 
degrees. While this offset may have initially been viewed 
as a strike against a DM annihilation interpretation of 
the line, our work demonstrates as a proof-of-principle 
that just such an offset should perhaps be expected. 

In Fig. [Tl] we show the DM annihilation luminosity 
"surface density" , J p^ M d£, calculated from the central 
(2 kpc) 3 region in Eris. The view is through the disk 
plane (the vertical axis corresponds to the disk normal) 
and the coordinate system has been rotated such that 
the offset peak is viewed edge-on, i.e. maximizing the 
angular offset. As before, the 3D density grid was first 
smoothed with a Gaussian kernel of width a = e so f t . We 
then squared it, multiplied the cells by their volume to 
get a luminosity, summed them up along the Y-axis, and 
finally divided the resulting map by the area of each cell 
to get a surface densitjj^jin units of M Q pc -5 . This map 
does not take into account any luminosity produced by 
DM along the line of sight outside of the central cube. 

The contrast between the maximum of the map 
(870M|pc -5 ) and the center (780M Q pc~ 5 ) is only 
about 11% at z = 0. Furthermore there is a second 
sub-dominant peak to the left of the center with a lumi- 
nosity surface density of 830 M Q pc -5 , only 5% less than 
the global maximum. If anything the contrast at z = 
is a bit of an outlier, as the mean contrast between offset 
peak and GC over the past 4 Gyr is only ~ 5% with a 
3% standard deviation. The greatest contrast measured 
over all outputs is 15%. Such small c ontrasts may not 
be com patible with the measurement of |Su fc F inkbeiner 
( 2012[ ), but it is important to keep in mind that the in- 
ternal structure of the peak is certainly not resolved in 
the simulation. It is conceivable that the contrast would 
increase with higher resolution, especially if resonances 
with the stellar bar are responsible, since those would be 
spread out artificially at low resolution. 

6. CONCLUSIONS 

We have analyzed the distribution of DM in the central 
regions of the hydrodynamical galaxy formation simula- 
tion Eris, one of the highest resolution and most realistic 
simulations to date of the formation of a barred spiral 
galaxy like our own Milky Way. Surprisingly, we find 
that the peak of the DM density in Eris is typically off- 
set from its dynamical center by several hundred parsec. 
No such offset is observed in its DM-only twin simulation 
ErisDark, nor in the much higher resolution DM-only Via 
Lactea II and GHalo simulations. 

5 To further clarify, the sum of all 200 2 "pixels" multiplied by 
the total area (4 X 10 pc 2 ) is equal to the total luminosity emitted 
from the (2 kpc) 3 cube. 
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The DM offset in Eris begins to appear around z — 1.5 
and grows over a period of 2 Gyr to a stable value of 
(-Doff) = 340 pc (almost three gravitational softening 
lengths), with a dispersion of 50 pc. The onset and du- 
ration of the DM offset appears to be well correlated 
with the formation of a nearly constant DM density core. 
The distributions of p max and D ff over the past 4 Gyr 
are inconsistent with a statistical fluctuation. Neither 
is the density peak a gravitationally bound structure, 
which rules out an incompletely disrupted subhalo core 
as an explanation. The most likely explanation may be 
a density-wave-like excitation by the stellar bar, possibly 
related to the resonant m echanism proposed by |Wein-| 
berg & Katz ( 2002 2007 ) to explain the transformation 
of a central DM cusp into a core. Arguments in favor of 
this explanation are the fact that the DM offset appears 
preferentially near the disk plane, that it is aligned to 
~ 30 degrees with the orientation of the stellar bar, and 
that is shows a periodicity of ~ 70 Myr. 

A central DM offset is of part icularly interest in the 
context of the recent report by Su & Finkbeiner ( 2012[ ) 
of a highly significant detection of gamma-ray line emis- 
sion from a region ~ 1.5 degrees (~ 200 pc projected) 
away from the Galactic Center. At first impression such 
a large angular offset would seem to argue against a 
DM annihilation interpretation of this signal. Our work 
demonstrates that in fact just such an offset should per- 



haps be expected. We note, however, that at the current 
resolution of our numerical simulations the low contrast 
(5 — 15%) of the annihilation surface brightness between 
the offset peak and the GC may be too small to accom- 
modate a DM annihilation explanation. 

We conclude by acknowledging that properly resolving 
the effects of baryonic physics on the central DM distri- 
bution, in particular those involving resonant bar-halo in- 
teractions, requires much higher resolution than we have 
been able to afford so far in cosmological hydrodynam- 
ics simulations. Further studies at higher resolution and 
exploring different baryonic physics implementation are 
sorely needed. Of particular importance are clarifying 
the role of the star formation threshold parameter and 
supernovae feedback for the formation of a DM core, as 
well as the influence of the supermassive black hole at 
the GC on the DM distribution. 
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